The isotropic-to-nematic transition in confined fiquid crystals : an essentially 

non-universal phenomenon 
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Computer simulations are presented of the isotropic-to-nematic transition in a liquid crystal con- 
Hned between two parallel plates a distance H apart. The plates are neutral and do not impose 
any anchoring on the particles. Depending on the shape of the pair potential acting between the 
particles, we find that the transition either changes from first-order to continuous at a critical film 
thickness H = H^, or that the transition remains first-order irrespective of H. This demonstrates 
that the isotropic-to-nematic transition in confined geometry is not characterized by any univer- 
sality class, but rather that its fate is determined by microscopic details. The resulting capillary 
phase diagrams can thus assume two topologies: one where the isotropic and nematic branches of 
the binodal meet ai H — H^, and one where they remain separated. For values of H where the 
transition is strongly first-order the shift Ae of the transition temperature is in excellent agreement 
with the Kelvin equation. Not only is the relation Ae (xl/ H recovered but also the prefactor of the 
shift is in quantitative agreement with the independently measured bulk latent heat and interfacial 
tension. 

PACS numbers: 64.70.M-, 05.70.Jk 



I. INTRODUCTION 



It is generally accepted that the first-order isotropic- 
to-nematic (IN) transition in liquid crystals confined 
between two parallel plates becomes continuous when 
the distance H between the plates becomes small 

Indeed, many simulations are consistent with this 
picture (7Hll| and show that the first-order IN transi- 
tion terminates at a critical film thickness H^. Some of 
these studies have also provided evidence of a continu- 
ous transition taking place when H < Hy^. Note that as 
H the system becomes effectively two-dimensional 
(2D). More recently, a (mathematically rigorous) proof 
appeared, showing that first-order IN transitions in 2D 
are also possible |12l . [Tsj . Inspired by this proof, com- 
puter simulations of liquid crystals in 2D were performed, 
which indeed uncovered strong first-order IN transitions 
too [3, [iHl- Hence, the IN transition in confinement 
can be continuous, as well as first-order. Finally, there 
is the scenario of no transition occurring at all in thin 
films [l^ [l3|, not even a continuous transition of the 
Kosterlitz-Thouless (KT) type |18]. Regarding experi- 
ments on confined liquid crystals, it has proved difficult 
to resolve continuous IN transitions in thin films [2, 19]. 
Pronounced coexistence between isotropic and nematic 
domains is typically observed p^| - |2l| . which suggests 
that a transition does take place and that it is first-order. 

The qualitatively different manifestations of the IN 
transition in confinement (continuous, first-order, ab- 
sence) rule out any universality class for this transition. 
What remains of the IN transition in thin films is de- 
termined by microscopic detail. The only regime where 
some "agreement" may be obtained is in the bulk 3D 
limit 7J — > oo. Here, one usually observes a first-order IN 
transition, with long-range order in the nematic phase. 



The transition thus breaks the rotational symmetry of 
the isotropic phase. At the mean-field level, this implies 
that the transition must be first-order |22|. We empha- 



size that fluctuations can change this result: even in 3D 
bulk, a genuine continuous IN transition is also possible 
(2^ [23} . However, most bulk experiments yield a first- 
order IN transition, and so the mean-field approximation 
appears to be valid in this regime. As the film thickness 
H decreases, fluctuations become increasingly important, 
and we expect three scenarios to unfold. In the first and 
most commonly accepted scenario, the IN transition be- 
comes continuous when the film thickness drops below 
a critical thickness H^. In addition, confinement is ex- 
pected to destroy long-range order in the nematic phase, 
due to the Mermin- Wagner theorem j^^j. Instead, quasi- 
long-range order may result, where the orientational cor- 
relations decay as a power law with distance. In the sec- 
ond (lesser known) scenario, the IN transition remains 
first-order irrespective of the film thickness, i.e. all the 
way down to — > 0. In the third scenario, evidence for 
which was recently provided [l^ [13 j the IN transition 
vanishes completely in the thin-film limit. 

Given the three scenarios for the IN transition in con- 
finement, all of which are qualitatively different, it is of 
fundamental interest to establish which "microscopic de- 
tail" is responsible for the scenario that ultimately oc- 
curs. The aim of this paper is to identify one possible 
mechanism, using computer simulations of a generalized 
Lebwohl-Lasher (LL) model. As it turns out, the gener- 
alized LL model is capable to reproduce all three scenar- 
ios, by tuning just a single parameter in the Hamiltonian. 
The effect of this parameter is to make the pair interac- 
tion "sharp and narrow" , meaning that particles interact 
when aligned but are otherwise rather indifferent to each 
other. Depending on this parameter, the crossover with 
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decreasing film thickness from first-order to continuous 
behavior can be ehminated completely, and the IN tran- 
sition remains first-order irrespective of H. 

The outline of this paper is as follows. We first in- 
troduce the generalized LL model and describe the sim- 
ulation method. Next, we present new simulation data 
showing one example where the IN transition becomes 
continuous below a critical film thickness H^, and a sec- 
ond example where the transition remains first-order ir- 
respective of the film thickness. We do not consider the 
scenario where the transition vanishes below as this 
has recently been done elsewhere [l^, [13] • A stringent 
test of the Kelvin equation, describing the shift of the 
transition temperature as a function of film thickness is 
also included. We end with a discussion and summary in 
Section HVl 



II. MODEL AND SIMULATION METHOD 

We consider a lattice model similar in spirit to the 
LL model [1^. To each site i of a 3D lattice, a 3D unit 
vector di (spin) is attached, which interacts with its near- 
est neighbors via 

E^-eJ2\d,-d,f, (1) 

with exponent p and coupling constant e. In this work we 
absorb a factor of l/ksT in the coupling constant, with 
ks the Boltzmann constant and T the temperature, and 
so e plays the role of inverse temperature. The lattice 
is a L X L X H rectangular box, with periodic boundary 
conditions in the lateral L directions but not in the H 
direction. The parameter H thus plays the role of the 
film thickness; the minimum thickness that can be stud- 
ied in this way equals H = 1, corresponding to a single 
lattice layer. This setup is identical to the slab geometry 
used in earlier simulations of the confined LL model Ji,. 
Note that spins at the walls have a lower number of near- 
est neighbors, but that the walls are otherwise neutral, 
i.e. we do not impose any anchoring conditions. 

In the original LL model the exponent of Eq. ([1]) equals 
Pll = 2. In the bulk limit H ^ oo, a (weak) first-order 
IN transition is observed ^Et^ at « 1.34 In 
the thin-film limit H — 1 recent simulations indicate the 
absence of any transition when p = 2 [l^, [l3| . In this 
work we consider p > 2. This modification is expected 
to enhance first-order phase transitions 13111 . which may 
then even survive the limit H ^ 1 [3, Il5j |. In line with 
previous work [23,[H,|3l|, we analyze Eq.([T]) in terms of 
the histogram 

P{E,S) = P{E,S\H,L,e), (2) 

defined as the probability to observe a system with energy 
E and nematic order parameter S' in a sample of thick- 
ness H , lateral extension L and at inverse temperature e. 



The distribution is obtained by computer simulations us- 
ing Wang-Landau [3l,[33| and transition matrix [s^ sam- 
pling; additional details pertaining to the present model 
are provided in Ref. [H. The nematic order parameter S 
is defined in the usual way as the maximum eigenvalue 
of the oricntational tensor 

1 ^ 

QaP = -^^^{SdiadjS - Sal3) , (3) 
i=l 

with dia the a component {a ~ x, y, z) of the orientation 
di of the spin at site i, the sum over all N = HL^ lattice 
sites and Sap the Kronecker delta. In a perfectly aligned 
sample it holds that S" = 1, whereas an isotropic sam- 
ple yields — > in the thermodynamic limit (hence, S 
defined in this way is an intensive quantity). A final in- 
gredient of this work is the use of finite-size scaling (FSS); 
needed because we seek thermodynamic limit properties. 
The thermodynamic limit of a film of thickness H is de- 
fined by extending the lateral extension L — > oo. In the 
bulk thermodynamic limit, both H and L are taken to 
infinity. 

III. RESULTS 

Depending on the exponent p in Eq. ([T]) , we expect the 
first-order IN transition either to terminate at a critical 
film thickness or to remain first-order irrespective of 
H . The case p = 2, i.e. the original LL model, is an 
example of the former scenario. In the bulk limit one 
obtains a first-order transition (25l - l29j , which terminates 
when the film thickness equals ^8—16 lattice lay- 
ers Q. In the 2D limit H ~ 1 no phase transition is 
observed for p = 2 |16|, ll7| . We emphasize that the lat- 
ter finding is not without some controversy, as previous 
other numerical studies of this system concluded that a 
phase transition does take place, namely a continuous 
transition of the KT type (see discussion in Ref. [TgI ) . 

A. crossover scenario 

We now consider Eq.(IT|) using a larger exponent, p = 8, 
to demonstrate that also a continuous IN transition is 
possible in thin films. To determine the order of the 
transition we use two FSS methods: the first was ini- 
tially given by Lee and Kosterlitz [36] and is based on 
the energy distribution P{E), defined as the probability 
to observe a system with energy E 

P{E) = j J S{E - E')P{E',S')dE'dS', 

with P{E, S) the joint distribution of Eq.©. 

At a first-order transition P{E) becomes bimodal, see 
Fig.[T]for an example, where the logarithm of the distribu- 
tion is shown. For finite L the bimodal structure persists 
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FIG. 1: Logarithm of P using p = 20 in Eq.lH) with H = I 
and 1/ = 25. The value of e has been chosen to give peaks of 
equal height. The free energy barrier, labeled AJ", is given 
as the difference between the peak maxima straddling the 
minimum. The distance labeled Ap corresponds to the latent 
heat density. The distribution is plotted as a function of the 
negative energy density, such that the left peak corresponds to 
the isotropic phase and the right peak to the nematic phase. 
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FIG. 2: Evidence of the crossover scenario, whereby the IN 
transition ceases to be first-order below a critical film thick- 
ness; the results in this plot refer to p = 8 in Eq.([l|. The 
free-energy barrier AF versus the lateral film extension L is 
plotted for several values of the film thickness H in (a). For 
large H the barriers increase linearly with L, consistent with a 
first-order transition; for smaller H the barrier vanishes with 
increasing L. The maximum value of the specific heat versus 
L is plotted in (b), again for several H. 



over a range of e values. As L increases the range be- 
comes smaller and in the thermodynamic limit L — ^ oo 
there is only one e where P{E) is bimodal, then featuring 
two 5-peaks. Hence, for finite L there is some freedom in 
choosing e and in Fig. [T] we have tuned e such that the 
peaks are of equal height. 

At a first-order transition the peak height, /S.F in 
lnP(i?), see the vertical arrow in Fig. [TJ corresponds to 
the free energy cost of interface formation [s^j ■ We there- 
fore expect /S.F oc L'^~^, with L the lateral extension of 
the film and d = 2 (recall that films are effectively two- 
dimensional). To determine the order of the transition, 
Lee and Kosterlitz [s^ proposed to measure AF versus 
L, which should yield a linear increase for a film. Re- 
sults for several values of H are shown in Fig. [JJa) . The 
data clearly indicate that the crossover scenario is taking 
place: for H ^ Q AF increases linearly with L, consistent 
with a first-order transition. In contrast, for H = 3 AF 
vanishes for large L, implying the absence of a first-order 
transition. 

To obtain the crossover thickness more accurately 
we use a second FSS method, based on the specific heat 



C={{E')~{Ef)/N, 



(4) 



with N — HL^ the number of lattice sites (volume). For 
given L and H a graph of C versus e reveals a maximum; 
the value of the maximum defines CL,ma.x{H). At a first- 
order transition the maximum scales with the volume of 
the system, that is CL,max{H) oc N [s^. In a film of fixed 
thickness H this implies CL,ma.x{H) oc L" with aist = 2. 
The result is shown in Fig. [Ifb) for several values of the 
film thickness. For = 6 a fit yields a = 2.00, con- 
firming that the transition is first-order. For H = 4 we 
obtain a = 1.74, indicating that a first-order transition is 
absent. Hence, we conclude that the crossover thickness 
ffx = 5. Precisely at a fit yields a = 1.94, which is 
still very close to the first-order value. Presumably for 
H — 5 the IN transition is weakly first-order. 

For H > iJx, i-e. where the transition is distinctly first- 
order, there is two-phase coexistence at the transition 
inverse temperature. It seems natural to characterize 
the phases with the nematic order parameter 5*. This 
approach is somewhat dangerous as confinement could 
destroy long-range nematic order in the thermodynamic 
limit: limi^oo 5* = irrespective of e. For H = 1 this 
follows rigorously from the Mermin- Wagner theorem ^24j] . 
The practical problem, affecting both simulations and ex- 
periments [1^ , is that the decay of S with L may be very 
slow. In fact, finite samples at low temperature typically 
reveal substantial order, even when the Mermin- Wagner 
theorem applies (39| . The present simulations are no ex- 
ception. Shown in Fig.[31[a) is S versus e in the bulk limit 

— > oo for several system sizes L (the bulk simulations 
were performed on a 3D cube of edge L with periodic 
boundaries in all directions). A first-order IN transition 
taking place at e w 1.52 [32], where S jumps to a finite 
value, is clearly seen. More importantly, for e above the 
transition, S becomes independent of system size, at least 
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FIG. 4: Capillary phase diagram of Eq. using p = 8. Shown 
is the variation of the coexisting phase densities piso{H) and 
Pnem(-ff) with the inverse film thickness 1/H. The critical in- 
verse thickness is at l/H^ ~ 0.2, above which the transition 
is no longer first-order and hence the two branches termi- 
nate. In the region between both branches coexistence be- 
tween isotropic and nematic phases is observed. 
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FIG. 3: Variation of the nematic order parameter S versus 
inverse temperature e using p = 8 in Eq.([I| for several val- 
ues of the film thickness H. In (a) we show the bulk result 
H — >■ oo, whereas (b) and (c) were obtained in films of finite 
thickness H. Note that for (a) and (b) the IN transition is 
first-order while it has become continuous in (c). 



on the scale of the graph; the latter is consistent with the 
formation of long-range nematic order, as expected in 3D. 
In Fig. IHb) we show the corresponding result for a film 
of thickness H — 10, which is still above the crossover 
thickness, and so the transition remains first-order. The 
behavior is similar to the bulk case, in the sense that S 
"jumps" at the transition, and for large e it appears to 



saturate at a finite value independent of the lateral film 
extension L. Hence, Fig. EJb) provides no evidence of 
S decaying to zero in the thermodynamic limit L — > oo, 
but rather that the film supports long-range nematic or- 
der. If S eventually does decay to zero, it is clear that 
huge system sizes, beyond the reach of any foreseeable 
simulation, are required to observe it. 

To avoid these subtleties, we characterize the coexist- 
ing isotropic and nematic phases in the film with their 
energy densities piso{H) and Pnom{H) respectively. These 
are simply the peak positions in the energy distribution 
of Fig. [1] Recall that the latent heat of the transition 
equals Cl{H) = Pnem{H) — piso{H), where the subscript 
is a reminder of finite-size effects in the lateral film ex- 
tension. The latent heat is related to the specific heat 
maximum (ssj 

£l{H) = ^4CL.nUH)/N (5) 

and the extrapolation to L oo is performed assuming 
that Coo {H)—Cl{H) cx The average energy density 

^^^^^ ^ pUH)+^PnUH) ^l_J ^p(^)^^ 

obtained at the specific heat maximum is extrapolated 
analogously: Pao{H) — pl{H) oc Once Coo{H) and 

Poo (H) have been determined, the coexisting energy den- 
sities follow. The latter may then be plotted in a capillary 
phase diagram, see Fig.|4l where the coexistence densities 
versus inverse film thickness l/H are shown. Since the 
transition ceases to be first-order at the critical thickness 
i/x the isotropic and nematic branches terminate. 
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FIG. 5: Cumulant analysis of Eq.([T]) using p = 8. Shown is 
Ui versus e using several values of the lateral film extension 
L, for (a) H = 1, (h) H ^ 2, and (c) H ^ 4. The value 
of e at the cumulant intersection yields the transition inverse 
temperature eoo{H) of the thermodynamic limit L ^ oo. 



We now consider H < H-^. For H ^ 1 and p = 2 in 
Eq. (IlJ , recent results [Tfil, flT'] indicate the absence of any 
phase transition (not even a continuous transition of the 
KT type). Part of the evidence is based on the failure of 
the Binder cumulant to intersect. At a continuous phase 
transition the ratio Ui = (S'^)/(5')^ becomes indepen- 
dent of system size [13, Slj , where S is the nematic order 
parameter. In simulations, this can be used to locate a 



TABLE I: Phase transition properties for H < H^, for the 
continuous IN transition. Listed is the transition inverse tem- 
perature eoa{H), along with the exponents j3 and 7, versus the 
film thickness H . The results refer to p = 8 in Eq.([T]). 
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continuous transition, by plotting Ui versus e for several 
system sizes L. At the transition inverse temperature 
too{H) of the film in the thermodynamic limit the curves 
for different lateral extensions L are expected to inter- 
sect. While for iJ = 1 and p = 2 no intersections are 
found [13], the result for p = 8 is radically different, see 
Fig. [S^a). Shown is J7i versus e using H = 1 for several 
values of L. The curves clearly intersect and so we con- 
clude that a continuous phase transition is taking place. 
This result strikingly illustrates the non-universality of 
the IN transition: whether a transition occurs for H — I 
is determined by the exponent p in Eq.([T]), i.e. a micro- 
scopic detail! Using p = 8 we have verified that contin- 
uous transitions exist for all values of the film thickness 
H < H^. The results for = 2 and H = 4, where H is 
approaching H^, are shown in Fig. [51^b) and (c), both of 
which reveal cumulant intersections. 

The fact that the cumulants intersect is a consequence 
of hyperscaling. At the transition inverse temperature 
(-oo{H) the order parameter decays (5) oc L^^ ^ while the 
susceptibility x = N{{S'^) - (Sf) diverges x oc L^f . The 
exponents /? and 7 are connected via the hyperscaling 
relation 



7 + 2/3 = d, 



(6) 



with spatial dimension d — 2 for a film. This relation im- 
plies that the order parameter and its root-mean-square 
deviation scale oc with the same exponent x. Conse- 
quently, appropriately constructed cumulant ratios, such 
as [/i, become independent of L whenever hyperscaling 
holds. By tuning the inverse temperature e we have de- 
termined (H) in our simulations by requiring that the 
scaling of (S) and x with L conforms to hyperscaling, 
i.e. we numerically solved Eq.®. A solution to Eq.® 
for each H < could indeed be found; the resulting es- 
timates of (-ff), as well as the exponents /? and 7, are 
listed in Table HI As expected, eoo{H) in Table |T] is close 
to the cumulant intersections of Fig. [Sj the discrepancy 
being less than 0.1 %. Note also that eoo{H) increases 
with decreasing H. The latter is consistent with the gen- 
eral tendency of confinement to lower phase transition 
temperatures. 

For H = 1 the system has become 2D and the expo- 
nents refiect "pure" values, free from any crossover ef- 
fects. Note that the exponents for H — 1 deviate signif- 
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FIG. 6: Variation of the susceptibility x with inverse temper- 
ature e for several values of the lateral film extension L, using 
film thicknesses {a) H = 1 and (b) H — 2. Both of these 
values H < and so the IN transition is continuous. Note 
the logarithmic vertical scale! The data were obtained using 
p = 8 in Eq.([l]). 



FIG. 7: Variation of the specific heat C with inverse temper- 
ature e for several values of the lateral film extension L, using 
film thicknesses (a) // = 1 and (b) H — 2. Both of these 
values H < and so the IN transition is continuous. The 
data were obtained using p = 8 in Eq.([T]). 



icantly from the XY values $xy = 1/8 and jxy = 7/4 
[4^ . strongly suggesting a different universality class. For 
H > 1, the trend is that /3 — > 0, while 7 — 2. Our 
interpretation is that, for 1 < H < H^, one observes 
crossover scaling behavior Q, governed by two compet- 
ing fixed points: one being the first-order transition at 
H = and the other being the continuous transition 
at = 1. The exponents for 1 < H < are there- 
fore "effective exponents" , with values between those of 
the H — 1 system, and the "first-order" values /3ist = 
and 7ist = d = 2 [38]. Note that effective exponents do 
not convey any fundamental information: if we were able 
to simulate arbitrarily large L values arbitrarily close to 
the transition inverse temperature, the same exponents 
as for the H = 1 system would be found. 

The important result to take from this analysis is that 
for p = 8 in Eq.([T]) and H < a continuous phase 
transition is found; by enforcing hyperscaling the transi- 
tion inverse temperature eoo{H) can be quite accurately 



obtained. 

We now consider how the nematic order parameter 5' 
depends on e and L; a typical result is shown in Fig.[3Jc) 
where H — 2 was used. We note that S increases with 
e and that the slope dS/de reaches a maximum close to 
^oo{H). In contrast to the first-order transitions observed 
for 7J > 7?x, S does not saturate at high e but decreases 
steadily with increasing L; this behavior is typical for all 
H < H^. Our simulation data thus suggest the absence 
of long-range nematic order in the thermodynamic limit 
when H < H^. This rules out a conventional critical 
point, since then the order parameter grows as a power 
law S (X , t > 0, implying S* > in the nematic phase, 
with distance from the transition 

t = e-eooiH) (7) 

and l3 the critical exponent of the order parameter. It is 
most likely, therefore, that the continuous transition we 
observe is a topological transition of the KT type [11] . 
Consistent with the KT scenario is our previous re- 
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suit of the order parameter decaying (S) oc L~^, and 
the susceptibihty diverging % oc L^, whilst obeying hy- 
perscaling. For completeness, we provide in Fig. [5] some 
raw simulation data for the susceptibility. Clearly visible 
is that X versus e reveals a maximum, becoming more 
pronounced for increasing L. In principle, the inverse 
temperature CL.xiH) where the susceptibility reaches its 
maximum, in a film of thickness H and lateral extension 
L, can be extrapolated using 



ln(L/c)iA 



(8) 



with non-universal constants b and c, and where the ex- 
ponent characterizes the exponential divergence of the 
correlation length ^ oc exp(5i'^) for t < 0, with t given 
by Eq.jT]). For the XY model it holds that vxy = 1/2, 
but since we did not recover XY exponents in Table H] 
the application of Eq.® requires that i/ be fitted also, 
implying a 4-parameter fit. We found that such a fit- 
ting procedure was numerically difficult to perform, and 
hence we did not determine eoo{H) in this manner. 

Finally, we note that also the specific heat, defined 
in Eq.Q, is consistent with the KT scenario. Plotted in 
Fig.[7]is the variation of C with e for several i, using two 
values of the film thickness. In both cases a maximum 
is revealed, but for _ff = 1 it grows only weakly with L. 
This is consistent with a negative specific heat exponent, 
implying that C remains finite in the thermodynamic 
limit, which agrees with the KT scenario. For H — 2 we 
observe that C already grows quite profoundly with L. 
We again attribute this to the crossover to a first-order 
transition where, ultimately, the specific heat maximum 
should scale oc L", with aist = 2, see also Fig.[2fb). 



B. first-order transitions 

We now consider the IN transition using p = 20 in 
Eq.(IT]). In this case, the transition is strongly first-order, 
even in the thin- film limit. The application of the Lee- 
Kosterlitz scaling method for H = 1 is shown in Fig.|51[a), 
where the linear increase of the barrier AF with L is 
clearly visible. The scaling of the specific heat maximum 
also confirms a first-order transition, see Fig.[8jb), show- 
ing the expected quadratic dependence of CL,raa^{H) on 
L. Since increasing the film thickness makes the transi- 
tion more strongly first-order, it is clear that for p — 20 
no crossover can occur. In the capillary phase diagram, 
see Fig. |9l the isotropic and nematic branches of the co- 
existing energy densities do not terminate, but continue 
all the way to H ^ 1. 



C. Kelvin equation 

Finally, we study the variation of the inverse transi- 
tion temperature too {H) with the film thickness for those 
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FIG. 8: Scaling analysis of Eq.lT]) using p = 20 and film thick- 
ness H = 1. In (a) we show the variation of the barrier AF 
versus L, while in (b) the specific heat maximum CL,max(^^) 
versus L is shown. Both these results indicate a first-order 
phase transition, even though the system is purely 2D. The 
dashed line in (a) is the result of a linear fit through the ori- 
gin. The curve in (b) is a fit to the form Cl, max (^^) oc L"; we 
obtain a. « 1.98, which is very close to fiist = 2 of a first-order 
phase transition in 2D. 




FIG. 9: Capillary phase diagram of Eq.([T]) using p = 20. 
Shown is the variation of the coexisting energy densities 
piao{H) and PnamiH) versus the inverse film thickness 1/H. 
In this case no crossover occurs and the IN transition re- 
mains first-order irrespective of _ff. The isotropic and nematic 
branches of the binodal therefore do not terminate, but con- 
tinue all the way to iif = 1. 
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TABLE II: Dependence of the transition inverse temperature 
^oo{H) on the film thickness H, for selected values of H where 
the IN transition is first-order. Results are shown for expo- 
nents p = 8 and p = 20 in Eq.([T]). The variation of eoo{H) 
with H should follow the Kelvin equation, see Eq.©. The 
bottom three lines list the bulk {H — oo) transition inverse 
temperature eoo the bulk latent heat density £oo , and the bulk 
interfacial tension , which are required in order to compare 
to the Kelvin equation. 



p = 8 H 




p = 20 


H 




5 


1.614 




1 


2.769 


6 


1.593 




2 


2.175 


7 


1.578 




4 


1.962 


8 


1.568 




8 


1.874 


10 


1.555 




10 


1.858 


15 


1.540 




30 


1.821 


30 


1.528 








50 


1.525 








100 


1.522 








Eoc 


1.521 




too 


1.806 




0.909 






1.727 


7oo 


0.06 




7oo 


0.30 



cases where the IN transition is first-order. We expect 
^oo{H) to fit to the Kelvin equation Q as 

Ae = l-e^/eoo(ff)- 



(9) 



where 700 is the bulk (iJ -> 00) interfacial tension, too 
the bulk IN transition inverse temperature and Coo the 
bulk latent heat density. In deriving this equation com- 
plete wetting is assumed All quantities that appear 
in Eq.® can, in principle, be extracted from finite-size 
simulation data with relative ease. For example, too{H) 
at a first-order transition can be obtained from ei^fc(iJ); 
the latter is defined as the inverse temperature where 
the ratio of the peak areas in the energy distribution 
P{E) equals k. For an optimal value k = fcopt, which 
can be found using trial-and-error, the L-dependence in 
(-L,k{H) becomes negligible and €oo{H) can be accurately 
obtained [s^l- The resulting estimates of the transition 
inverse temperatures, for both p = 8 and p = 20, are pro- 
vided in Table |lll using only values of the film thickness 
where the transition is first-order. 

Similar to previously, bulk iJ — > 00 results are obtained 
using L X L X L systems with periodic boundaries in all 
directions. The bulk latent heat density Coo is obtained 
from the specific heat maximum using Eq. ([5]) and is once 
again extrapolated to L 00, where now N — L^. The 
resulting estimate of £00 is also listed in Table [Til To 
obtain the bulk interfacial tension 700 we use the method 
of Binder 



37 



Simulating a large and stretched LxLxD 
system D > L with periodic boundaries in all directions, 
the logarithm of the energy distribution P{E) reveals a 
pronounced fiat region between the peaks, see Fig. [TUl 




FIG. 10: Plots of InP as obtained in completely periodic 
simulation boxes of size 10 x 10 x 30 (solid lines) and 10 x 10 x 60 
(dashed lines) for (a) p = 20 and (b) p = 8. The height 
of the peaks AF is related to the interfacial tension 70c via 
Eq. (|10|) . The distance Ap between the peaks is a measure 
of the latent heat density £00- In these plots e was tuned to 
yield an approximately horizontal region between the peaks. 



The flat region indicates that the isotropic and nematic 
phase coexist with only small interactions between the 
two interfaces. Hence, the average peak height AF is 
related to the bulk interfacial tension 



7o 



hm 7i, 7i=AF/(2i2), 



(10) 



yielding an elegant method of obtaining 700. Provided 
L is large enough, the result should not depend on the 
elongation D, but inspection of Fig. [TUl reveals this is not 
quite true, especially for p = 8. This could indicate that 
some interaction between the interfaces remains, or that 
L was not large enough. In any case, using the largest 
available system size, we obtain ~ 0.05 for p = 8 and 
7oo ~ 0.29 for p = 20 (in units of /c^T per lattice spacing 
squared). Alternatively, 7^, can be measured in a cubic 
periodic system of size L and extrapolation to L — > cxd 
using 



7l = 7oo + ci In i/L^ + C2/L^ 



(11) 



with constants ci, can be attempted [SJ]. When using 
this procedure we obtain slightly higher values of the in- 
terfacial tension, namely 700 ~ 0.08 and 700 ~ 0.31 for 
p = 8 and p = 20 respectively. Hence, for p — 20 the es- 
timates for 7oo agree reasonably well, whereas for p = 8 
some discrepancy remains. In Table |ll] the average of 
both estimates is provided. 

We now have all quantities needed to put the Kelvin 
equation to the test, see Eq.®. Shown in Fig. [Tl]is Ae 
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0.20 




FIG. 11: Test of the Kelvin equation. Plotted is the inverse 
temperature shift Ae of Eq.© versus the inverse film thick- 
ness 1/H , using exponents p = 8 (a) and p = 2Q (b) in Eq.([T|). 



versus 1/H, for both p = 8 and p — 20, using only val- 
ues of H where the transition is first-order. Provided the 
Kelvin equation holds, the resulting plots should be lin- 
ear. For p = 8 this is clearly not the case; only in the 
limit 1/ H — > 0, i.e. where the transition is strongly first- 
order, is agreement observed. In contrast, using p — 2Q 
the Kelvin equation holds for all values of the film thick- 
ness, including H — 1. The slope a of the lines in Fig.fTT] 
can be obtained from a fit; following Eq.® it is expected 
that a = 27oo//^oo, allowing for a stringent quantitative 
test. For p = 20 we obtain by fitting a w 0.34, which is in 
excellent agreement with I'^ooj C.oo ~ 0.35 calculated us- 
ing the independent estimates of Table [TTl For p = 8 the 
fit yields a ~ 0.14, where only the largest three values of 
H were used. Once again, this is in excellent agreement 
with 27oo/i3oo ~ 0.13 obtained from Table HIl 



IV. DISCUSSION AND SUMMARY 

In this paper we have provided new results regarding 
the IN transition in liquid crystals confined between neu- 



tral walls. The main conclusion to be taken from this 
work is that a single universal scenario describing the 
nature of this transition as function of the film thickness 
R does not exist. Using a generalized version of the LL 
model, we have explicitly demonstrated that the first- 
order IN transition can terminate at a critical thickness 
i/x, below which it becomes continuous, or that it can 
stay first-order irrespective of H . The scenario that takes 
place is determined by a single parameter in the Hamil- 
tonian, namely p in Eq.([T]), which sets the "sharpness" of 
the pair interaction. When the transition is sufficiently 
strongly first-order excellent agreement with the Kelvin 
equation is also obtained. In particular, we not only ob- 
serve the 1 jH shift of the transition inverse temperature 
but also the prefactor of the shift is in quantitative agree- 
ment with the independently measured bulk latent heat 
and interfacial tension. However, when the IN transition 
is only weakly first-order clear deviations appear and the 
Kelvin equation significantly underestimates the inverse 
temperature shift, see Fig. fTTTa). 

The two different manifestations of the confined IN 
transition presented in this work yield two distinct phase 
diagram topologies: one where the isotropic and nematic 
branches of the binodal terminate at the critical thickness 
and one where they continue irrespective ol H. It is of 
some interest to compare the resulting phase diagrams to 
other works. The topology of the p = 8 phase diagram, 
see Fig. |4l is commonly encountered in confined colloidal 
rods and plates [ol-fTTi [isj . To facilitate the comparison, 
the energy density in Fig. 3] should be interpreted as the 
analogue of the particle density in colloidal systems. In 
agreement with Fig. |4l the first-order IN transition in 
colloidal systems also terminates at a critical thickness 
[ol-fTTI, l43j . It is also interesting to see that the nematic 
branch of the binodal in Fig. 2] shows rather extreme out- 
ward curvature as the bulk limit is approached. Colloidal 
platelets reveal similar behavior, albeit that here the ef- 
fect appears in the isotropic branch [43] . In contrast with 
colloidal systems is the fact that Eq.([T]) with p = 8 in the 
bulk limit yields a first-order transition that is too strong. 
Defining the relative strength of the transition as 



(12) 



Pnem ~r Piso 



we obtain r « 0.38 for Eq.([T|) with p = 8, while Onsager's 
exact solution (IJj for infinitely slender rods yields r « 
0.12. This discrepancy can be fixed by usin g a lower p 
in Eq.(IT]). For instance, p = 5 gives r « 0.15 |32l |. which 
is much closer to Onsager's result. Note that p = 5 still 
exceeds the original LL value p = 2. Indeed, it has been 
pointed out that the original LL model yields a bulk IN 
transition that is too weakly first-order compared to what 
is observed in fluids of rods [45| . 

The second phase diagram topology, where the bin- 
odal branches do not terminate in thin films, is obtained 
for p = 20 in Eq.(IT]), see Fig. [HI The resulting phase 
diagram is of fundamental importance, since it clearly 
demonstrates that first-order IN transitions in thin films 
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are also possible and that the crossover to a continuous 
transition need not necessarily take place. It is inter- 
esting that experiments so far have not produced clear 
evidence of a continuous IN transition in thin films [2l.[T9l- 
[2l| . This is consistent with a phase diagram topology as 
shown in Fig. HI However, it is obvious that the model 
of Eq.([T|) with p = 20 does not capture the bulk limit 
correctly, since the bulk IN transition ought to be weak, 
whereas p = 20 yields a very strong first-order transition. 
Clearly, some features are still lacking in Eq.(IT]), for in- 
stance a coupling between the orientational and spatial 
degrees of freedom of the particles, as well as anchor- 
ing effects at the walls. Investigations which incorporate 
these effects are possible directions for future work. 

Finally, using p = 8 and H < our results show 
that a genuine continuous IN transition can also take 
place. Since long-range nematic order is not observed a 
transition of the KT type [T^ is the most likely scenario. 



This result is interesting because using p = 2 one finds 
that Eq.(IT|) is without any kind of phase transition in 
the thin-film limit [l^, [l3]. Hence, the nature of the 
IN transition in thin films is ultimately determined by 
microscopic details. This means that a single universality 
class for the IN transition cannot exist. Depending on 
the details of the interaction, there can be both first- 
order and continuous transitions as well as no transition 
occurring at all. 
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